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Abstract 

New parameter-sets of the M3Y-type semi-realistic nucleon-nucleon interaction are obtained, by 
taking into account the whole contribution of the interaction to the pairing. The interactions are 
applicable to nuclear structure studies via mean- field calculations and their extensions. Imple- 
menting self-consistent mean-field calculations for the spherical nuclei, we confirm that the results 
do not change significantly from the previous ones, in which the parameters were determined by 
partly discarding influence of the density-dependent repulsion on the pairing. 
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I. INTRODUCTION 



Mean-field (MF) theories provide us with a basic tool to study nuclear structure from 
nucleonic degrees of freedom. However, it is yet difficult to reproduce structure of medium- 
to heavy-mass nuclei with the bare nucleon-nucleon (NN) interaction to good accuracy, 
and most of the MF calculations have been performed with phenomenologically determined 
effective interactions (or energy density functionals). 

The Michigan 3-range Yukawa (M3Y) interaction [l| was obtained by fitting the Yukawa 
functions to Brueckner's G-matrix. As the crucial role of the chiral symmetry breaking has 
been recognized, it is desirable to take into account the pionic effects in the NN interaction. 
The M3Y interaction contains the central part of the one-pion exchange potential (OPEP) 
and the tensor part that is corrected due to the medium effects as well as to the p-meson 
exchange. The author has developed semi-realistic effective NN interactions 0,0], by mod- 
ifying the M3Y interaction so as to reproduce the saturation properties and the is splitting 
while keeping the central part of the OPEP. The tensor channels have been pointed out 
to be significant in the Z- or iV- dependence of the shell structure We have obtained 
parameter-sets in which the tensor channels are maintained (e.g. the set M3Y-P5), and one 
in which the tensor channels are dropped (M3Y-P4) [3] . By implementing MF calculations, 
it has been shown that semi-realistic interactions could give different shell structure from 
the widely used Skyrme and Gogny interactions, and that of the semi-realistic interaction 
including the tensor channels is often favorable if compared to the experimental data. Sig- 
nificance of the tensor channels in the magnetic excitations has also been confirmed [5] via 
calculations in the random-phase approximation. 

The parameter-sets developed in Ref. jsl were obtained by fitting them to pairing proper- 
ties, as well as to data of doubly magic nuclei. However, in the calculations of finite nuclei in 
Ref. [3J , influence of the density-dependent contact force on the pairing was not fully taken 
into account. Because we find that this influence is not negligibly small, we propose modified 
parameter-sets in this report, which will be useful in future studies of nuclear structure. The 
new parameter-sets are refitted to the same quantities as the old ones, as will be described 
in the subsequent section. 

II. M3Y-TYPE INTERACTION 

We consider the effective NN interaction that has the following form, 
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TABLE I. Parameters of M3Y-type interactions. 
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where = — r,, = | r^- 1 with i and j representing the indices of nucleons, and 
p(r) denotes the nucleon density. P$e, Pte, Pso an d -Pro denote the projection operators 
on the singlet-even (SE), triplet-even (TE), singlet-odd (SO) and triplet-odd (TO) two- 
particle states. The density-dependent contact force uf? has been added to reproduce 
the saturation |. The Yukawa fuuetiou ff>(,r) = e^/^r is assumed for X = C, 
LS and TN. The modified parameter-sets M3Y-P4' and P5' are presented in Table [J in 
comparison with the previous sets M3Y-P4 and P5. The non-central channels and v^J 
are not changed; v$ S) and v { ™ ] of M3Y-P4' (M3Y-P5') are identical to those of M3Y-P4 
(M3Y-P5) U The M3Y-P5' interaction has vjj^ determined from the G-matrix jsj, while 
v ij^ — is assumed in M3Y-P4'. Not displayed in Table [Q the longest range (n = 3) part of 
v\j is kept to be f OPEP (the central channels of the OPEP). The other strength parameters 
are fitted to the measured binding energies of 16 and 208 Pb and to the even-odd mass 
differences of the Sn isotopes, in the Hartree-Fock (HF) or the Hartree-Fock-Bogolyubov 
(HFB) approximation. Compared to the old sets M3Y-P4 and P5, in M3Y-P4' and P5' 
the SE channel is more attractive below the saturation density, and its influence on the 
saturation and the symmetry energies is compensated by the other central channels. 
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TABLE II. Nuclear matter properties at the saturation point. 
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III. NUCLEAR MATTER PROPERTIES 



Energy of the nuclear matter is a function of p = J^arPra, Vs — J2ar vPral P, Vt — 
J2ctt T Pra/ P and rj st = J2ctt crT Pra/ P hi the HF approximation, with p Ta (r = p,n and a =t, 4) 
stands for densities depending on the spin and the isospin. Formulas to calculate the nuclear 
matter energy and its derivatives for given p Ta have been derived in Ref. 

The energy per nucleon £ = E/A has its minimum £ Q , at the saturation density po 
(equivalently, kpo). The incompressibility and the volume symmetry energy are defined by 

d 2 £ 



K = V 



dp 2 



ld 2 £ 



(2) 



2 dr# 

Here | indicates evaluation at the saturation point. The curvatures with respect to r) s and 
rjst are denoted by a s and a st . The fc-mass Mq is defined by a derivative of the s.p. energy. 
Density-dependence of the symmetry energy is typically represented by 
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These quantities calculated from the new semi- realistic interactions are tabulated in Table HT1 



For comparison, the values obtained by the D1S [7( and DIM [8] parameter-set of the Gogny 
interaction are also displayed. We also present the Landau-Migdal (LM) parameters in 
Table [TTTl which are calculated by using the formulas given in Ref. (2]. See Ref. Q also for 
definition of the LM parameters. It has been known that g' ~ 1 0], which is reproduced by 
the semi-realistic M3Y-type interactions primarily due to VQp E p. 

For the symmetric nuclear matter, M3Y-P4' and P5' give £{p) very close to each other, 
and also to those of M3Y-P4 and P5 j3j. £(p) in the neutron matter is shown in Fig. [TJ As 
in M3Y-P4 and P5, the p-dependence of the SE channel in M3Y-P4' and P5' gives sizable 
repulsion at high p, by which £{p) is qualitatively similar to the result of Ref. [10|, unlike 
the D1S case. 
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TABLE III. Landau-Migdal parameters at the saturation point. 
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IV. APPLICATION TO FINITE NUCLEI 

The HF and HFB calculations are implemented in finite nuclei by using the Gaussian 



expansion method [llHl4j and adopting the Hamiltonian H = Hp? + Vc — H c m , where H^, 
Vc and H c _ m , represent the effective nuclear Hamiltonian, the Coulomb interaction and the 
center-of-mass Hamiltonian. The exchange term of Vc is treated exactly. Both the one- and 
the two-body terms of H c , m , are subtracted before iteration. 

The binding energies and rms matter radii of several doubly magic nuclei are calculated 
in the spherical HF approximation. The results of the new semi-realistic interactions are 
compared with those of D1S and DIM, as well as with the experimental data, in Table ITVl 
The new interactions reproduce the binding energies and the rms matter radii with similar 
accuracy to M3Y-P4 and P5. The deviation in the binding energies is ^ lOMeV. Though 
slightly worse than D1S and DIM, we do not take this deviation seriously, until correlations 
due to the residual interaction are taken into account. 

In M3Y-P4' and P5', some of the parameters in the SE channel are fitted to the even- 
odd mass differences A^ ass (iV) = E(Z,N) - \[E{Z,N + 1) + E(Z,N- 1)] (N = odd) of 
the 66 < iV < 80 Sn isotopes, in the spherical HFB approximation with the £ < 7 basis 
truncation. Contribution of vlf D ^ to the pairing is fully taken into account. The equal filling 
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FIG. 1. Energy per nucleoli £ = E/A in the neutron matter. Blue solid, brown dotted, green dot- 
dashed and red solid lines represent the results with the D1S, DIM, M3Y-P4' and P5' interactions, 
respectively. Circles are the result of Ref. [lol ]. 



TABLE IV. Binding energies and rms matter radii of several doubly magic nuclei. Experimental 



data are taken from Refs. 15l4l8|. 
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FIG. 2. Even-odd mass difference in the Sn isotopes, A^= 5 s °(iV). The results of D1S, M3Y-P4' 
and M3Y-P5' are shown by blue open circles, green open squares and red diamonds, respectively. 
Experimental values, presented by black crosses, are taken from Ref. jl8j ]. 



approximation 19] is applied to the odd-mass nuclei. See arguments in Ref. [3j on influence 
of the I truncation, the particle-number conservation and the non-spherical mean fields. 
A^g°(iV) calculated by M3Y-P4' and P5' are compared with the experimental data and 
with those obtained by D1S in Fig. El Though not shown, DIM gives A^ S g°(iV) quite close 
to that of D1S except at iV = 92. The mass difference at N ~ 64 and 90 reflects interaction- 
dependence of the shell structure. In the M3Y-P4' and P5' interactions the N = 64 subshell 
effect is stronger than in D1S. At iV ~ 90 the M3Y-type interactions yield larger mass 
difference than D1S, because nlf 7 / 2 and n2p 3 / 2 well mix due to the pairing. 

We depict difference between the HF and the HFB energies, which represents the pair 
correlation, for the Ca nuclei in Fig. [3] and for the Ni isotopes in Fig. HI These results 
are qualitatively similar to those shown in Ref. [jj. However, we dare to repeat several 
interesting points. The vanishing difference between the HF and HFB energies is normally 
connected to the shell closure. It is found that 52 Ca is nearly a doubly-magic nucleus with 
any of the interactions, as is consistent with the experiments (2^]. With M3Y-P5', 68 Ni is 
almost doubly magic in harmony with experimental data 21[, although 60 Ca is not. This 
suggests significant Z-dependence of the shell structure around N = 40, whose origin will be 
argued elsewhere. The hindrance of the pair excitation at 86 Ni suggests magic or submagic 
nature of N = 58 due to the gap between n2s\/2 and nld^/2- 

Prediction of the neutron drip line depends on effective interactions to a certain degree. 
Compared to the results with M3Y-P4 and P5 [sj], more neutron-rich Ca and Ni nuclei are 
bound with M3Y-P4' and P5'. The heaviest bound Ca (Ni) nucleus has iV = 50 (66) in both 
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FIG. 3. Difference between the HF and HFB energies for the Ca isotopes (N = even), obtained 
from D1S, M3Y-P4' and P5'. See Fig. [2] for conventions. Dotted lines are drawn to guide eyes. 




FIG. 4. Difference between the HF and HFB energies for the Ni isotopes (JV = even). See Fig. [2] 
for conventions. 



S 



of the M3Y-P4' and P5' results. 
V. SUMMARY 

We have obtained new parameter-sets of the M3Y-type semi-realistic effective interac- 
tion to describe structure of nuclei, by taking into account the whole contribution of the 
interaction to the pairing in the HFB calculations. Two new parameter-sets have been pre- 
sented; one keeping the tensor force of the M3Y-Paris interaction (M3Y-P5'), and the other 
discarding the tensor force (M3Y-P4'). The OPEP part in the central force is maintained. 
Basic characters of the interactions are examined in the infinite nuclear matter and in several 
doubly magic nuclei. The pairing properties are checked by the even-odd mass difference of 
the Sn isotopes in the HFB approximation. Pair energies and location of the neutron drip 
line are investigated for the Ca and Ni isotopes. The results do not change significantly from 
the previous ones, in which influence of the density-dependent contact force on the pairing 
was partly discarded. 
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